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Phase diagram of an Ising-spin Kondo lattice model on a triangular lattice near 1/3-filling is 
investigated by Monte Carlo simulation. We identify a partially disordered phase with coexistence 
of magnetic order and paramagnetic moments, which was unstable in two-dimensional Ising models 
with localized spins only. The partial disorder emerges in the competing regime between a two- 
sublattice stripe phase and three-sublattice ferrimagnetic phase, at finite temperatures above an 
electronic phase separation. The peculiar magnetic structure accompanies a charge order and de- 
velops a gap in the electronic structure. The results manifest a crucial role of the nonperturbative 
interplay between spin and charge degrees of freedom in stabilizing the partial disorder. 

PACS numbers: 75.30.Kz,75.10.-b,75.40.Mg 
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The antiferromagnetic (AF) Ising model on a trian- 
gular lattice is one of the most fundamental models for 
geometrically frustrated systems. When the interaction 
is restricted to the nearest-neighbor (n.n.) pairs, frus- 
tration in each triangle prevents the system from long- 
range ordering (LRO) down to zero temperature (T), and 
the ground state has extensive degeneracy and associated 
residual entropy [lHl] . The degenerate ground state is ex- 
tremely sensitive to perturbations; e.g., further-neighbor 
interactions can lift the degeneracy and induce a variety 
of LRO. 

An intriguing state emergent from the degenerate 
ground-state manifold is a partially disordered (PD) 
state. The PD state is peculiar coexistence of magnet- 
ically ordered moments and thermally-fluctuating para- 
magnetic moments. Such possibility was first discussed 
in the presence of next n.n. ferromagnetic (FM) inter- 
action if. A mean- field study reported the existence of 
a three-sublattice PD phase with an AF ordering on the 
honeycomb subnetwork and paramagnetic moments at 
the remaining sites [Fig. QJd)]. Although such PD state 
was indeed observed in several Co compounds @ , sub- 
sequent Monte Carlo (MC) simulations indicated that, in 
purely two-dimensional systems, the PD appears at most 
as a quasi-LRO and that the transition from the high- 
T para mag netic phase is of Kosterlitz-Thouless (KT) 
type |7H13|. The results suggest that, for establish- 
ing a PD LRO, it is indispensable to incorporate addi- 
tional elements, such as a three-dimensional interlaycr 
coupling [13 1 . 

In this Letter, we explore the possibility of PD LRO 
in two dimensions when taking account of the coupling 
to itinerant electrons. Our study is partly motivated by 
the recent discovery of PD in metallic compounds with 
quasi- two-dimensional structure [14j. The interplay be- 
tween localized moments and itinerant electrons plays a 
crucial role in the following points. First of all, the ki- 
netic motion of electrons induces effective interactions 
known as the Ruderman-Kittel-Kasuya-Yosida (RKKY) 
mechanism 15l- 17|. The long-ranged and oscillating na- 



ture of interaction drive keen competition between dif- 
ferent magnetic states. Furthermore, the change of mag- 
netic states affects the electronic state in a self-consistent 
manner through the spin-charge coupling; the system 
can gain the energy by forming some particular elec- 
tronic state associated with magnetic ordering. We note 
that similar PD states were recently discussed as a con- 
sequence of quantum spin fluctuation or hybirdization 
between itinerant and localized electrons jlil . Il9j ]. In 
the present study, however, we will show that a mini- 
mal model, in which neither the quantum spin fluctua- 
tion nor hybridization operates, does exhibits a PD LRO. 
Our results will provide deeper understanding of PD phe- 
nomena robustly observed in a broad range of materials, 
not only in the insulating localized spin systems but also 
the systems including itinerant electrons. 

We consider a single-band Kondo lattice model on a 
triangular lattice with localized Ising spin moments. The 
Hamiltonian is given by 



(i) 



The first term represents hopping of itinerant electrons, 
where Ci a (cJ CT ) is the annihilation (creation) operator of 
an itinerant electron with spin a =f , \. at ith site, and t 
is the transfer integral. The sum is taken over n.n. 
sites on the triangular lattice. The second term is the 
onsite interaction between localized spins and itinerant 
electrons, where of represents the z-component of itiner- 
ant electron spin and Si = ±1 denotes the localized Ising 
spin at ith site; J is the coupling constant (the sign of J 
does not matter in the present model). Hereafter, we take 
t = 1 as the unit of energy, the lattice constant a = 1, 
and the Boltzmann constant fee = 1. To investigate ther- 
modynamic properties of this model, we adopted a MC 
simulation which is widely used for similar models poj ]. 
The calculations were conducted up to the system size 
N = 18 2 under the periodic boundary conditions. The 
thermal averages were calculated for typically 4300-9800 
MC steps after 1700-5000 MC steps for thermalization. 



To describe the PD and fcrrimagnetic (FR) states, here 
we use the pseudo-spin defined for each three-site unit 
cell, 
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and its summation M = J2 m 3S m /AT, where m is the 
index for the three-site unit cells, and (i,j,k) denote 
the three sites in the mth unit cell belonging to the 
sublattices (A,B,C), respectively Then, the three- 

sublattice PD state [Fig. |TJd)] is characterized by a fi- 
nite M = (M x ,M y ,M z ) parallel to (yjij2, 1 1^/2,0), 
(0, V% 0), or their threefold symmetric directions around 
the z-axis. On the other hand, the three-sublattice FR 
state [Fig. [He)] is characterized by a finite M along 
{y/2/3, \/2, 1/V3), (2^/2/3, 0, -1/Vt), or their threefold 
symmetric directions around the z-axis. Hence, the two 
states are distinguished by the azimuth of M in the 
xy-planc (as well as M z ). To parametrize the situa- 
tion, we introduce ip = M^ y cos64>M, where 4>m is the 
azimuth of M in the xy plane and M xy = 3M xy /8 
(M xy = M x + My). The parameter ip has a negative 
value and ip for the PD ordering, while it be- 

comes positive and ip — > 1 for the FR ordering; ip = for 
both paramagnetic and KT phases in the bulk limit. We 
calculate ip, M and corresponding susceptibility, and spin 
entropy S (defined later) to map out the phase diagram. 

Figure [lj a) shows the phase diagram around the elec- 
tron density n — J2 ilT ( c l a c ia) /N = 1/3 at J = 2 ob- 
tained by MC calculations. There are four dominant or- 
dered phases — the stripe, PD, FR, and KT phases, in 
addition to an electronic phase separation (PS). In the 
relatively low density region n < 0.30, the stripe order 
with period two [Fig. QJc)] develops in the low T region. 
This is an interesting state, in which the magnetic order 
breaks the sixfold rotational lattice symmetry; Reflect- 
ing the lower symmetry, the electronic transport becomes 
spatially anisotropic. On the other hand, in the higher 
density region n > 1/3, the system exhibits the three- 
sublattice FR order [Fig. [2(e)] at low T. In this den- 
sity region, we find another phase transition at a higher 
T than the onset of FR order. As described below, we 
identify the intermediate phase as a KT-like state sim- 
ilar to the KT state discussed in the Ising models 0- 
11] . In the intermediate region 0.30 < n < 0.34, the 
three-sublattice PD phase [Fig.QJd)] emerges at finite T, 
whereas it is taken over by PS (or FR) at low T. In the 
finite-T region, the numerical data indicate that the PD 
state retains LRO, in sharp contrast to the KT region, 
as discussed below. 

Figure [2] shows typical T dependences of physical quan- 
tities in each density region. Figures EJa) and [2Kb) 
present the results in the stripe region at n = 0.27. 
Figure [U(a) shows the order parameter for stripe or- 
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FIG. 1. (color online). Phase diagram of the model |l} (a) 
at J = 2 while varying n and (b) at n — 1/3 while varying 
J. The symbols shows phase boundaries for the four phases: 
stripe phase, partially-disordered phase (PD), KT-like phase 
(KT), and fcrrimagnetic phase (FR). PS represents a phase 
separation. The lines are guides for the eyes. The colored 
strips at T — show the ground states obtained by comparing 
the energy of stripe and FR orders. The boundary between 
PD and PS in (b) is lower than T = 0.02, and difficult to 
determine by MC calculations. The schematic pictures of the 
magnetic orders are given for (c) stripe, (d) PD, and (e) FR 
phases. The circles in (d) show paramagnetic sites. 



der, M str = E q . {S(q s * tr )/AT} 2 ] 1/2 , where 5(q) is 
the spin structure factor of Ising spins and q* tr = 
(7r, 0), (±^7r, -^7r) are the characteristic wave vectors of 
the stripe order; Figure d^b) represents the correspond- 
ing susceptibility Xstr- A phase transition to the stripe 
phase is signaled by a rapid increase of M s t r and corre- 
sponding peak of Xstr j 

the peak T of Xstr, T c (stl) , for each 
system size is plotted in the phase diagram in Fig. HJa). 

The results in the PD region for the pseudo-spin mo- 
ments and the corresponding susceptibilities at n = 1/3 
are shown in Figs.[^c) and[U(d), respectively. M xy shows 
a rapid increase around Tc™' = 0.130(4) and approaches 
y/2 at low T, whereas |M Z | does not grow. Correspond- 
ingly, Xxy shows a divergent peak at Tc PD K The low 
temperature phase is identified as the PD phase by the 
azimuth parameter ip shown in Fig.[3Ia); ip becomes neg- 
ative below Tc PD ^ and approaches — 1|, as expected for 
the PD instability. At low T, the PD phase becomes un- 
stable and is taken over by PS, which is determined by 
the jump of n as a function of the chemical potential. 

Figures[Ue) and[^f ) show the results in the FR region 
at n = 0.38. The data indicate two successive transitions 



3 









( a ) - 




'\ \-Mstr 










□ 7V=12xl2 




o iV = 12xl8 




*7V=18xl8 









1.0 
0.8 
0.6 
0.4 
0.2 



0.0 

0.03 0.04 



1.2 
1.0 
0.8 
0.6 
0.4 
0.2 



0.05 

T 





(c) - 






- « co [*30]*||». 


Myi/^i - 


_ i ,V-12il2 % 




o JV=12xl8 




_ ',¥=18x18 










cassis; 



0.06 0.07 0.03 

180 
150 
120 
90 
60 
30 




0.07 







" *\ X z [x3] 






fa H v' 












$< 





0.0 

0.00 0.04 0.08 0.12 0.16 0.20 0.00 0.04 0.08 0.12 0.16 0.20 



1.2 
1.0 
0.8 
0.6 
0.4 
0.2 





fa, M xv [/V2]^ e ) 












" □ A" =12x12* 




_oA<=12xl8 




ijV = 18xl8 



0.0 

0.04 0.08 



0.12 

T 



0.16 0.20 0.04 




0.20 



FIG. 2. (color online). MC results for (a) M str , (b) Xstr, 
(c),(e) M xy , \M Z \, and nco, and (d),(f) % xy and Xz for the 
system sizes N = 12 x 12, 12 x 18, 18 x 18. The data are at 
(a),(b) n = 0.27, (c),(d) n = 1/3, and (e),(f) n = 0.38. 



indicated by the peaks in \ xy and \ z at different T. The 
peak of Xz with corresponding increase of \M Z | signals the 
phase transition to the FR phase at T^™' = 0.098(4). 
On the other hand, at a higher Tkt = 0.146(4), only 
M xy changes rapidly, and correspondingly, Xxy shows a 
peak. M xy , however, shows a noticeable system-size de- 
pendence even below Tkt, in contrast with the results 
below Tj PD ^ in Fig. [2jc). Similar behavior was observed 
in the KT ordering in Ising spin systems ficlllll ]. Further- 
more, ip does not show any clear anomaly at Tkt, while 
it shows a sharp rise around Tj FR ^ , as shown in Fig. [3Jb) . 
The value of ip extrapolated to large N converges to zero 
in the intermediate T range. This indicates that there is 
no sixfold symmetry breaking in M xy at Tkt, as seen in 
the KT phase in the Ising spin models [l(| ■ The extrap- 
olated ip becomes finite below Tj FR ' and approaches 1 as 
expected for the FR ordering. Hence, we consider that 
the higher-T transition is of KT type; the system exhibits 
two successive transitions from the paramagnetic phase 
to the KT phase, and the KT phase to the low-T FR 
phase. 

The PD phase in the intermediate n exhibits charac- 
teristic features reflecting the presence of the paramag- 
netic moments. Figures Etc) andE^d) show the results of 
effective spin entropy S [21[ at n = 1/3 and 0.38, respec- 
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FIG. 3. (color online). MC results for (a),(b) the tj) param- 
eter and (c),(d) spin entropy S and its derivative dS/dT for 
the system sizes TV = 12 x 12, 12 x 18, 18 x 18. The data are at 
(a),(c) n = 1/3 and (b),(d) n = 0.38. S are calculated from 
numerical integration of DS/dT, assuming S(T = 1) = log2. 



tively. The data for dS/dT are also shown, which give 
the specific heat associated with the Ising spins divided 
by T. For n = 1/3 [Fig. [3{c)], S shows a steep decrease 
around T ~ Tj P , corresponding to the phase transi- 
tion to the PD state. However, it does not decrease to 
zero below Tc and stays at a value close to I log 2, the 
residual entropy expected for 1 /3 disordered moments in 
the ideal three-sublattice PD state [22[ . In contrast, S 
in the FR region [Fig. [3]f d)] shows little change at Tkt, 



whereas a sharp decrease is seen at T, 



(FR). 



S becomes 



smaller than ^ log 2 and approaches zero while decreasing 
T, as expected for the FR LRO. The contrasting behavior 
depending on the presence and absence of the paramag- 
netic moments is also seen in the susceptibility Xz shown 
in Figs.^d) andUff). 

We also note that the PD state accompanies concomi- 
tant charge modulation. In Fig.[2jc), we show the result 
for the charge order parameter nco = {^(cico)/-^} 1 ^ 2 
at q^Q = (— 27r/3, 27r/\/3) [-AT(q) is the charge struc- 
ture factor for the itinerant electrons]. The result shows 

a rapid increase of nco below T<; PD ', indicating that 
a three-sublattice charge order develops with PD (local 
charge density becomes poor at the paramagnetic sites). 

Performing the analyses above in a wide range of pa- 
rameters, we conclude that the system exhibits the three- 
sublattice PD LRO in the region near n = 1/3 at finite 
T, while it turns into a KT-like quasi-LRO for larger n, 
as summarized in the phase diagram in Fig. [IJa). We 
also extend the analyses while changing J and confirm 
that the PD state near n = 1/3 remains stable in the 
wide range of J. For instance, Fig. QJb) shows the phase 
diagram at n = 1/3 while varying J; The PD state ap- 
pears for J > 1.0, while the FR as well as KT state is 
stabilized for smaller J. 
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FIG. 4. (color online). MC result for DOS of itinerant 
electrons at n = 1/3 for N = 18 x 18. The Fermi level is 
set at e = 0. The solid curve shows the integrated DOS at 
T = 0.04. 



Let us discuss our results in comparison with the previ- 
ous ones for the Ising spin models. As shown in Fig.QJa), 
the PD state appears in the phase competing region be- 
tween the stripe and FR phases. Although similar phase 
competition occurs in the Ising spin models while chang- 
ing the second-neighbor interaction in the presence of the 
dominant AF interaction for n.n. pairs, PD does not ap- 
pear in the form of LRO; the system shows, at most, a 
quasi-LRO with a KT transition 0-E3- The PD LR0 
in our model, however, emerges in the distinct param- 
eter regime from the KT region; it dominantly appears 
above the PS region in the phase competing regime. PS 
is an electronic origin associated with a finite jump of n 
at the magnetic transition, and does not appear in the 
Ising spin models. Hence, our results indicate that the 
electronic degree of freedom opens a new window where 
the PD can be stabilized [23j |. 

In addition, Fig.[IJb) shows that the PD phase appears 
only above a finite J and is not stable down to the small J 
limit. This indicates that the RKKY interactions derived 
by the second-order perturbation in terms of J [l5- 17 
are not sufficient to explain the appearance of PD state. 
The result implies that the nonperturbative effect beyond 
the RKKY interactions takes a crucial role in stabilizing 
the PD LRO " 



24] 



Related to this point, we find that the PD state is 
stabilized by the change of the electronic structure. Fig- 
ureUJshows the density of states (DOS) at n = 1/3 while 
varying T. The result shows the development of an en- 
ergy gap at the Fermi level while decreasing T. The gap 
structure is always observed in the PD region, whereas 
such feature is not seen in the KT region. We note that a 
'mean-field' PD state with assuming the three-sublattice 
order of (+to, — m, 0) can open a charge gap for a suffi- 
cient J. These results clearly indicate that the PD LRO 



is stabilized by an energy gain from the reconstruction of 
electronic state with gap opening. 

To summarize, we have provided the convincing nu- 
merical results on the partial disorder in the Isin-spin 
Kondo lattice model on a triangular lattice. Interestingly, 
our partially-disordered state emerges at finite tempera- 
tures dominantly above an electronic phase separation 
inherent to the itinerant model. Furthermore, we found 
that the partial disorder develops with a charge gap in the 
electronic structure as well as charge ordering. These re- 
sults indicate that the nonperturbative interplay between 
spin and charge degrees of freedom is crucial for stabiliz- 
ing the partial disorder. Experimentally, recent studies 
on quasi-two-dimensional metallic oxides reported an in- 
dication of partial disorder. Our results will promote 
understanding of the interesting physics related to the 
peculiar coexistence of magnetic order and paramagnetic 
moments in itinerant electron systems. 

The authors are grateful to H. Kawamura, M. Mat- 
suda, and H. Yoshida for fruitful discussions. The au- 
thors also thank S. Hayami and T. Misawa for help- 
ful comments. Part of the calculations were performed 
on the Supercomputer Center, Insitute for Solid State 
Physics, University of Tokyo. H.I. is supported by Grant- 
in-Aid for JSPS Fellows. This research was supported 
by KAKENHI (No.19052008, 21340090, and 22540372), 
Global COE Program "the Physical Sciences Frontier", 
the Strategic Programs for Innovative Research (SPIRE), 
MEXT, and the Computational Materials Science Initia- 
tive (CMSI), Japan. 



[1] 

[3 
[5 



m 

[7] 

[«] 

[9 
[10 

[11 

[12 
[13 
[14 



G. H. Wannier, Phys. Rev. 79, 357 (1950). 

R. M. F. Houtappel, Physica (Amsterdam) 16, 425 
(1950). 

K. Husimi and I. Syozi, Prog. Theor. Phys. 5, 177 (1950). 
M. Mekata, J. Phys. Soc. Jpn. 42, 76 (1977). 
T. Kohmoto, T. Goto, S. Maegawa, N. Fujiwara, Y. 
Fukuda, M. Kunitomo, and M. Mekata, Phys. Rev. B 
57, 2936 (1998). 

S. Niitaka, K. Yoshimura, K. Kosuge, M. Nishi, and K. 
Kakurai, Phys. Rev. Lett. 87, 177202 (2001). 
K. Wada, T. Tsukada, and T. Ishikawa, J. Phys. Soc. 
Jpn. 51, 1331 (1982). 

S. Fujiki, K. Shutoh, Y. Abe, and S. Katsura, J. Phys. 

Soc. Jpn. 52, 1531 (1983). 

D. P. Landau, Phys. Rev. B 27, 5604 (1983). 

H. Takayama, K. Matsumoto, H. Kawahara, and K. 
Wada, J. Phys. Soc. Jpn. 52, 2888 (1983). 

S. Fujiki, K. Shutoh, S. Katsura, J. Phys. Soc. Jpn. 53, 
1371 (1984). 

T. Takagi and M. Mekata, J. Phys. Soc. Jpn. 64, 4609 
(1995). 

N. Todoroki and S. Miyashita, J. Phys. Soc. Jpn. 73, 412 
(2004). 

M. Matsuda, C. de la Cruz, H. Yoshida, M. Isobe, and 
R. S. Fishman, Phys. Rev. B 85, 144407 (2012). 



5 



[15] M. A. Ruderman and C. Kittel, Phys. Rev. 96, 99 (1954). 

[16] T. Kasuya, Prog. Theor. Phys. 16, 45 (1956). 

[17] K. Yosida, Phys. Rev. 106, 893 (1957). 

[18] Y. Motome, K. Nakamikawa, Y. Yamaji, and M. Uda- 

gawa, Phys. Rev. Lett. 105, 036403 (2010). 
[19] S. Hayami, M. Udagawa, and Y. Motome, J. Phys. Soc. 

Jpn. 80, 073704 (2011). 
[20] S. Yunoki, J. Hu, A. L. Malvezzi, A. Moreo, N. Furukawa, 

and E. Dagotto, Phys. Rev. Lett. 80, 845 (1998). 
[21] M. Udagawa, H. Ishizuka, and Y. Motome, Phys. Rev. 

Lett. 108, 066406 (2012). 
[22] The small deviations of ip and <S from their expected 



values are presumably ascribed to a spatial correlation 
between the disordered spins. 

[23] PS might be replaced by an inhomogeneous state by 
forming domains when taking account of additional el- 
ements such as impurities and Coulomb repulsions. Nev- 
ertheless the PD is expected to survive, because mag- 
netically, it is an Ising-type order robust against (non- 
magnetic) impurities, and electronically, it accompanies 
a charge order, which can be further stabilized by the 
long-range Coulomb repulsion. 

[24] We also conducted calculations by changing n, and found 
that the PD state always appears from a finite J. 



